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Abstract 

A quantum characteristic exponent may be denned, with the same 
operational meaning as the classical Lyapunov exponent when the lat- 
ter is expressed as a functional of densities. Existence conditions and 
supporting measure properties are discussed as well as the problems 
encountered in the numerical computation of the quantum exponents. 
Although an example of true quantum chaos may be exhibited, the 
taming effect of quantum mechanics on chaos is quite apparent in 
the computation of the quantum exponents. However, even when the 
exponents vanish, the functionals used for their definition may still 
provide a characterization of distinct complexity classes for quantum 
behavior. 
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1 Introduction. Classical and quantum char- 
acteristic exponents. 

A notion of quantum characteristic exponent has been introduced in Ref. [p]], 
which has the same physical meaning as the corresponding classical quan- 
tity (the Lyapunov exponent). The correspondence is established by first 
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rewriting the classical Lyapunov exponent as a functional of densities and 
then constructing the corresponding quantity in quantum mechanics. The 
construction is explained in detail in Ref.0, where the required functional 
spaces are identified and the infinite-dimensional measure theoretic frame- 
work is developed. Here we just recall the main definitions and emphasize 
some refinements concerning the support properties of the quantum charac- 
teristic exponents, which turn out to be relevant for the numerical computa- 
tions of Sect. 2. 

Expressed as a functional of admissible L 1 — densities, the classical Lya- 
punov exponent ispl 



A„ = lim — log 

n— >oo yi 



~ vi ~L* Ds * (/ My) y pnp{y) 



(i) 



where p is an initial condition density, P the Perron-Frobenius operator, x a 
generic phase-space coordinate, v a vector in the tangent space, \i the invari- 
ant measure and D$ x the Gateaux derivative along the generalized function 
8 X . The possibility to define Gateaux derivatives along generalized functions 
with point support and the need for a well-defined cr-additive measure in 
an infinite-dimensional functional space lead almost uniquely to the choice 
of the appropriate mathematical framework, that is, admissible densities are 
required to belong to a nuclear space. Being ergodic invariants, the Lya- 
punov exponents exist on the support of a measure. In the nuclear space 
framework, measures with support on generalized functions, which are in 
one-to-one correspondence with the usual measures in phase space, may be 
constructed by the Bochner-Minlos theoremfQ. 

To construct, in quantum mechanics, a quantity with the same opera- 
tional meaning as ([!]) let U n (n continuous or discrete) be the unitary opera- 
tor of quantum evolution acting on the Hilbert space H and X a self-adjoint 
operator in H belonging to a commuting set S. For definiteness X is assumed 
to have a continuous spectrum, to be for example a coordinate operator in 
configuration space. One considers, as in the classical case, the propagation 
of a perturbation di5 x , where by x we mean now a point in the spectrum of 
X. 

v'Dm (u n V, XU n ®) = 2Re A < £ U~ n XU n ^ > (2) 
1 v ' ox 1 
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For the proper definition of the right-hand side of (|2|) one requires ^ G E to 
be in a Gelfand triplet 

E* D if D £ 

By < 5^1 or < x | we denote a generalized eigenvector of X in E* . Notice also 
that U n , being an element of the infinite-dimensional unitary group, has a 
natural action both in E and -E"*[|J. One obtains then the following definition 

of quantum characteristic exponent 



\ V)X = lim sup — k>£ 

n— >oo Tl 
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dx l 
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The support properties of this quantum version of the Lyapunov expo- 
nent have to be carefully analyzed. In Eq.(|3]), ^ defines the state which, in 
quantum mechanics, plays the role of a (non-commutative) measure^. The 
quantum exponent may depend on the state, but the measure that, as in the 
classical case, provides its support is not the state but a measure in the space 
of the perturbations of the initial conditions, that is, in the space where the 
Gateaux derivative operates. In the classical case these two measures coin- 
cide, in the sense that to which invariant measure in phase-space corresponds 
an infinite-dimensional measure in the space of generalized functions @. In 
the quantum case, however, the two entities are different, the second one 
being the measure on the spectrum of X induced by the projection-valued 
spectral measure and the state, that is 

u(Ax) = J dP x m) (4) 

Ax 

A particular case where an infinite-dimensional measure-theoretical set- 
ting, similar to the classical one, may be used to define the quantum exponents! 
is when the quantum evolution is implemented by substitution operators 
in configuration space, as in some sectors of the configurational quantum 
cat H [|]] . However this formulation is not very useful in general and the 
state plus spectrum-measure framework seems to be the one that has general 
validity. In this framework the following existence theorem holds 

Theorem : Let X be a self-adjoint operator, E a test function space in a 
Gelfand triplet containing the generalized eigenvectors of X in its dual E* 
and * G E. Then if U n $ G E and X$ G E (V$ G E, Vn G Z) and the 
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following integrability condition is fulfilled 



dv{x) log 



\Rev l d Xi < x \U 


>| 









< M 
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V$ G E, the limit in Eq.@ exists as a L 1 (z/)— function, that is , the average 
quantum characteristic exponent is defined for any measurable set in the 
support of v. 



Proof: 

We write Eq.(OI) as 



i n-i rwa^ < x\u- n+k xu n - k <& > 

K(x) = lim sup — log — (6) 

«-»n n f Rev l d Xi < x\U~ n+k+l XU n - k - 1 ^ > 

Then from the integrability condition (|5|) the integral of the sequence in the 
right-hand side of Eq.(^) is bounded and the Bolzano- Weierstrass theorem 
insures the existence of the lim sup. Therefore X v (x) is well defined as a 
L l (v)— function. □ 

Notice that we really need the lim sup in the definition of the characteristic 
exponent because we have no natural [/-invariant measure in E to be able, for 
example, to apply Birkhoff or Kingman's theorem and prove limsup=liminf. 
Also the sense in which the measure v provides the support for the quantum 
characteristic exponent is different from the classical ergodic theorems. We 
have not proven pointwise existence of the exponent a. e. in the support a 
measure. What we have obtained here is the possibility to define an average 
quantum characteristic exponent for arbitrarily small ^—measurable sets. 

Other definitions of characteristic exponents in infinite-dimensional spaces 
have been proposed by several authors || [(J [|l(|] [|TTJ [|13|]. They charac- 



terize several aspects of the dynamics of linear and non-linear systems. The 
definition discussed here, proposed for the first time in seems however to 
be the one that is as close as possible to the spirit of the classical definition 
of Lyapunov exponent. 

Like the classical Lyapunov exponent the quantum analogue (|j) cannot 
in general be obtained analytically. There is however a non-trivial exam- 
ple where it can. This is the configurational quantum cat introduced by 
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Weigert [0] . The phase space of this model is T 2 x R 2 . A mapping similar 
to the classical cat operates as a quantum kick in the configuration space 
T 2 , and the rest of the Floquet operator is a free evolution. This system has 
the appealing features of actually corresponding to the physical motion of a 
charged particle on a torus acted by an impulsive electromagnetic field and, 
as show by Weigert , to be exactly solvable. 
The Floquet operator is 

U = U F U K (7) 
Up = exp[— i— p 2 }; Uk = exp[— -(x ■ V • p + p ■ V • x)\ (8) 



Up is a free evolution and Uk a kick that operates in a simple manner on 
momentum eigenstates and on (generalized) position eigenstates 

U K \p) = \m~ x p) (9) 

U K \x) = \Mx) (10) 

M being an hyperbolic matrix with integer entries and determinant equal to 
1. The momentum has discrete spectrum, p e (2irZ) 2 . 

To compute the quantum characteristic exponent (Eq.(|3|)), let the oper- 
ator X be 

X = sin(2vr/ ■ x) (11) 

I G Z 2 . This operator has the same set of generalized eigenvectors as the 
position operator x. To construct the measure v (Eq.(f|)) in the spectrum 
of the operator X we cannot use the energy eigenstates | Pa) because they 
are not normalized. However all one requires is invariance of the measure, 
and using the (normalizable) momentum eigenstates one such measure is 
obtained. 

u(A) = {p | / dx\x){x\p) (12) 

J A 

This invariant measure in this case happens to be simply the Lebesgue mea- 
sure in T 2 . 
Defining 

ln (x) = (x | U- n XU n | p) (13) 
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the result for the quantum characteristic exponent is 



A„ = limsup^^ \ log + Rev 1 ^ n {x) 



(14) 



= limsup n ^ 00 ilog + \v i (2irM n l) i {cos9 n (pJ,x) + cos0 n (p, -l,x)}\ 
with 

rp /n—1 n—l \ 

6 n (p, l,x) = -[ J2( M ~ k P) 2 + E (M k (2nl + M"») 2 + x • (2nM n l + p) 

2 \fc=0 k=0 ) 

For the lim sup the cosine term plays no role and finally 

A„ = lim sup - log v\M n l)i (15) 

n^oo TL 

The characteristic exponent is then determined from the eigenvalues of the 
hyperbolic matrix M and is the same everywhere in the support of the mea- 
sure v. If Hi (a*i > 1^2) are the eigenvalues of M, one obtains = log/xi 
for a generic vector v and A^ = log H2 iff ^ is orthogonal to the eigenvector 
associated to /ii- Hence, in this case, one obtains a positive quantum char- 
acteristic exponent whenever the corresponding classical Lyapunov exponent 
is also positive. 

This exact example will be used in Sect. 2 as a testing ground for the 
numerical algorithm and an illustration of the kind of precision problems and 
support properties to be expected when computing quantum characteristic 
exponents. 

In the numerical calculation of the quantum characteristic exponents two 
delicate points are identified. The first is that the calculation requires a high 
degree of precision, because, if the exponent is positive, the derivative of 
U~ n XU n ^(x) grows very rapidly with n. Therefore in the positive exponent 
case acceptable statistics is only obtained by taking average values over the 
configuration space. Second, if the situation is as in the classical case where 
different invariant measures coexist in phase space, the quantum exponent 
may depend on the state \P used to define the measure v in the spectrum of 
X. Therefore, in all rigor, one should first construct stationary states and 
then study the ^/—dependence of the quantum exponent. Such study has 
not yet been carried out and, in the calculations below, a flat wave function 
is used as the initial state. 
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2 Numerical computation of quantum char- 
acteristic exponents 

For kicked quantum systems corresponding to the Hamiltonian 

H = H Q + V{x) y £6{t-jT) (16) 

3 

the Floquet operator is 

U = e -iV{x) e -irH ^) (17) 

in units where Ti = 1. For the computation of the action of U on a wave 
function if)(x), a fast Fourier transform algorithm F and its inverse F _1 are 
used 

Uij)(x) = e- iV{x) F- l e- iTHo{lk) Fip(x) (18) 

In this way one obtains a uniform algorithm for any potential. In the con- 
figurational quantum cat, Eq.(H), the computation is similar with the multi- 
plicative kick e~* v ( x ' replaced by the substitution operator 

i/j(x) -> j/j(M -1 x) (19) 

The quantum characteristic exponent is obtained from the calculation of 

d x U~ n XU n ij{x) (20) 

in the limit of large n. The precision of the algorithm is checked by insuring 
that 

(y- n U n - l) il){x)\ <e (21) 
for a small e, and that the finite difference used to compute the derivative 



in (|20"D does not approach the maximum possible value allowed by the dis- 



cretization. 



2.1 The configurational quantum cat 

Here the configuration space is the 2-torus T 2 , the Floquet operator is the 
one of Eq.(|8p, and the kick is a substitution operator with matrix 

M={\\) (22) 
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Numerically we have computed the quantities 

- n XU n V) (x) \ 

(23) 




Re (lJ- n XU 




(*) 




Re 


(a) 





T2 



X being the operator in (|TT|). The initial wave function is = 1 and the 
average is taken over the whole of configuration space. It turns out that, in 
this case, the derivative in the numerator of (EBf) g rows so feist £it some points 
that one reaches, after a few iterations, the maximum difference (2 in this 
case) for the wave function at two nearby points in the discretization grid. 
When this happens the calculation cannot be reliably taken to higher n with 
that discretization. In the numerical calculation of the classical Lyapunov 
the computation becomes a local evaluation at each step by rescaling the 
transported tangent vector. Here, because of the linearity of matrix elements 
and quantum evolution, a similar procedure is not possible and one has to 
carefully control the growth of the quantities in (|23| ) . To improve statistics 
the average over configuration space has been taken. This can be safely done 
in this case because we know exactly the supporting measure (|12"D, but in 
general there will be no guarantee that the supporting spectral measure is 
uniform. In any case average quantities like (|23|) are exactly we expect to be 
able to compute reliably. 

Fig.l shows the evolution of - (D n — Dq) obtained with a discretization 
grid of 292681 points in the unit square, for two different directions v. The 
calculation was interrupted when the local finite differences reached one half 
of the maximum. The lines are fits to the points constrained to approach 
the same value at large n. The resulting numerical estimate for the largest 
quantum characteristic exponent is 0.95. The exact value obtained from (|15|) 
is 0.9624. 

2.2 Quantum kicked rotators 

The configuration space is the circle S 1 , 

V(x) = q cos(2vrx) (24) 
x G [0, 1) and for H the following two possibilities were explored 

H W = l^L 
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The operator X is 

X = sin(27re) (26) 
The quantity that is numerically computed is 




and, in all cases, one starts from a flat initial wave function. In both cases and 
for the very many values of q that were studied, this quantity seems either to 
stabilize or to have a very small rate of growth for large n. Fig.2, for example, 
shows the results for the Hq case with r = ^ and q — 5. The (numerical) 
conclusion is that the quantum characteristic exponent vanishes. This con- 
clusion does not seem to be a numerical artifact because the discretization 
grid for the fast Fourier transform has always been chosen sufficiently small 
to insure a small local finite difference for all iterations. For example in the 
example shown in Fig.2, the grid has 4096 points which keeps observed lo- 
cal differences below 0.1. Also the vanishing of the quantum characteristic 
exponent in quantum kicked rotators is not dependent on the phenomena of 
localization because also for quantum resonances it may exactly be shown to 
vanish [JJ. 

In Fig.2 (D n ) seems to tend to a constant at large n. In other cases very 

(2) 

slow rates of growth are observed. This is shown in Figs.3a,b for the Hq 
case with r = ^ and q — 11. 

3 Conclusions 

Both the classical Lyapunov exponent (HD and its quantum counterpart @, 
measure the exponential rate of separation of matrix elements of X when 
the density (or the wave function) suffers a S x perturbation, x being a point 
in the spectrum of X. The configurational quantum cat example shows 
that there are instances of true quantum chaos, in the sense of exponential 
growth of the matrix element separation. However, as the numerical study 
of the quantum kicked rotators seems to show, exponential growth may be 
rather exceptional in quantum mechanics. Furthermore the taming effect of 
quantum mechanics on exponential chaos goes deeper than the phenomenon 
of localization, because also for quantum resonances, where no localization 
is present, the quantum characteristic exponent vanishes. 
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Although distinct from one another, all known ways that now exist to 
approach the problem of quantum chaos, seem to agree in one point, namely 
that quantum mechanics has a definite taming effect on chaos. This is now 
probably the main issue in quantum chaos, not only from the theoretical 
point of view, but also in the context of quantum control. Even if quantum 
characteristic exponent, as defined in (|3|) might be zero in most quantum 
systems, the rate of growth index D n or its average (D n ) might still be 
useful as a characterization of quantum dynamics because, even if weaker 
than exponential, a growth of this quantity would still be an indication of 
sensitivity to initial conditions. In particular, as suggested by the numerical 
results, subexponential rates of growth might characterize distinct complexity 
classes of quantum evolution. 

4 Figure captions 

Fig.l - Calculation of ^ {D n — D ), Eq. (p3|) , in the configurational quantum 
cat for two orthogonal directions v and a fit constrained to the same limit at 
large n. 

Fig. 2 - (D n ), Eq . (f27D , for the quantum standard map at q = 5, r = 
Fig. 3 - (a) (D n ), Eq. (|27|) , for a kicked rotator with kinetic Hamiltonian 
Hq 2 ^ at q = 11, t = ^ ; (b) the same scaled by log(log(n + 1)). 
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